Vertex corrections in the dynamic structure factor in spin ladders 
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We combine the results of perturbative continuous unitary transformations with a mean-field 
calculation to determine the evolution of the single-mode, i.e., one-triplon, contribution to the 
dynamic structure factor of a two- leg 5 = 1/2 ladder on increasing temperature from zero to a finite 
value. The temperature dependence is induced by two effects: (i) no triplon can be excited on a 
rung where a thermally activated triplon is present; (ii) conditional excitation processes take place if 
a thermally activated triplon is present. Both effects diminish the one-triplon spectral weight upon 
heating. It is shown that the second effect is the dominant vertex correction in the calculation of 
the dynamic structure factor. The matrix elements describing the conditional triplon excitation in 
the two-leg Heisenberg ladder with additional four-spin ring exchange are calculated perturbatively 
up to order 9. The calculated results are compared to those of an inelastic neutron scattering 
experiment on the cuprate-ladder compound La4SrioCu2404i showing convincing agreement for 
established values of the exchange constants. 

PACS numbers: 75.40.Gb, 78.70.Nx, 75.10.Jm, 74.72. Cj 



I. INTRODUCTION 

Currently, the properties of strongly correlated systems 
such as quantum spin systems are much better under- 
stood in their respective ground states, i.e., at zero tem- 
perature, than at finite temperatures. Generically, the 
ground state is a pure state and the excitations may or 
may not have unconventional statistics. But since they 
occur only in small numbers at zero temperature their 
statistics is not exceedingly difficult to describe. 

At finite temperatures, however, the physics becomes 
richer and more complicated because there is a finite con- 
centration of excitations. Then their statistics and their 
mutual interactions play an important role. Sharp 5- 
peaks signalling infinitely long-lived excitations acquire 
a broadening induced by thermal fluctuations. Spectral 
weight is shifted due to variations of the matrix elements 
and due to variations of the available phase space. 

For one-dimensional quantum spin systems Essler and 
collaborators analytically investigated the changes of line 
shapes in the dynamic structure factor [TH3]. They con- 
sidered systems where the elementary excitations are 
hardcore bosons which are local, i.e., they have mini- 
mal spatial extension. The main finding is that the line 
shapes broaden upon heating as expected. The positions 
of the peaks shift depending upon the available phase 
space. Rather unexpectedly, the line shapes at finite tem- 
peratures are generically asymmetric, in agreement with 
earlier numerical evidence [4]. Recently, the experimen- 
tal verification of these findings has been achieved [5] . So 
the usual expectation of Lorentzian broadening does not 
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hold. This is due to the interplay of the harcore property 
and the available phase space. 

In the present article, we aim at the investigation of 
the complementary effect induced by vertex corrections. 
We focus on the changes in the spectral weight of the 
elementary modes due to heating. We do not study the 
broadening of the peaks corresponding to the elementary 
modes. 

The system we are investigating is a two-leg ladder 
with localized S — 1/2 at each vertex, see Ref. and 
references therein. The dominant couplings are the an- 
tiferromagnetic Heisenberg exchange couplings between 
nearest neighbours, see Fig. [I] The nearest neighbour ex- 
change couplings are of similar magnitude. The ground 
state has zero spin and the elementary excitations are 
of triplet character, i.e., they have S = 1. Thus we 
dub them triplons [7] in order to distinguish them from 
the elementary spin wave excitations, magnons, of sys- 
tems with magnetic long-range order. The centers of the 
triplons can be seen to sit on the rungs of the ladder. 
Further details about the system and its experimental 
realization will be given in the following sections. 

The temperature development of the spectral weight 
is induced by two effects. First, the triplons are hard- 
core bosons so that they can only be created on a rung 
if there is no other triplon on this rung before. Clearly, 
this effect depends strongly on temperature because the 
relative frequency of a successful excitation depends on 
the concentration of thermally activated triplons. The 
higher the temperature is the less probable is an excita- 
tion process. This effect, due to the hardcore nature, has 
been dealt with in the preceding investigations as well 

mm- 

Second, triplons can be excited if there are other 
triplons on rungs nearby. We call this process a con- 
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Figure 1: (Color online) Sketch of the two-leg spin ladder. 
Circles indicate spins with S — 1/2. Solid lines represent the 
two-spin exchange couplings on the rungs J± and on the legs 
J||. The four-spin cyclic exchange is denoted by J cyc . 



ditional excitation henceforth. Thereby, the matrix ele- 
ment for the excitation of a triplon depends on the con- 
centration of thermally activated triplons. It turns out 
that the sign of the conditional excitations is negative 
so that higher temperatures and concomitantly higher 
concentrations of triplons imply a lower spectral weight. 

The conditional excitations become more and more im- 
portant the more the triplons are extended. Hence the 
vertex corrections due to conditional excitations are sig- 
nificant for isotropic spin ladders where the couplings 
along rungs and legs of the ladder are of similar mag- 
nitude. In the limits investigated previously the vertex 
corrections did not play a significant role [TJ [51 [3] . 

We find that the conditional excitations dominate the 
temperature dependence of the spectral weight. Our the- 
oretical results agree very well with our experimental 
ones so that evidence for the relevance of these vertex 
corrections is provided. 

Technically, our theoretical calculations are based first 
on continuous unitary transformations [5J [S] which map 
the original model onto an effective model which con- 
serves the number of excitations [TUl - fT2"] . This is done in 
a perturbative fashion where the unperturbed Hamilto- 
nian is given by isolated rungs P3J Q3] . Alternative real- 
izations of degenerate perturbation theory are also based 
on systematic changes of the underlying basis [I5lfP7] . 

Second, we simplify the operators which describe exci- 
tation processes by a mean- field (MF) treatment. In this 
step we focus on a description in terms of single modes. 

The paper is set up as follows. In Sect. [TT] we present 
our model in detail. Subsequently, the theoretical ap- 
proach is presented in Sect. |III[ The theoretical results 
are shown in Sect. |TV] and compared to experimental re- 
sults in Sect. |V} Finally, the paper is concluded in Sect. 
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II. MODEL 



We study a two-leg 5=1/2 ladder as it is depicted 
in Fig. [TJ The model is characterized by Heisenberg ex- 
change couplings on the rungs J± and on the legs Jii of 
the ladder as well as by a four-spin cyclic exchange J cyc . 



The corresponding Hamiltonian in units of Jj_ reads 

Jii „ J c 

h 



H 



^CyC TT 

t -"eye 



la) 



(lb) 



where i denotes the rungs and r € {1,2} the upper and 
lower leg of the ladder. The expansion parameters em- 
ployed are 



^cyc 



:= Ju/Jj. 

■ — Jcyc I J A_ • 



(2a) 
(2b) 



The four-spin ring exchange is given by 

i 



(3) 



These four-spin terms constitute the most important cor- 
rection to the commonly considered nearest neighbor 
Heisenberg exchange. By now, this is established for 
planes [TBH2T] and ladders CM25]. The size of J cyc 
in cuprate ladders was found to be of the order 0.2-0.25 
J_i_ [26- 29J . Further corrections like the two-spin interac- 
tion over the diagonal are almost one order of magnitude 
smaller [53] so that they are neglected in the following. 

The physics of spin ladders is well understood from the 
limit of isolated rungs {x — 0,a; C y C =0}. In this limit 
the ground state is the product state of singlets and the 
excitations are completely local triplets. The elemen- 
tary excitations at finite values x and cc cyc are dressed 
triplets which are called triplons [7]. Triplon excitations 
are gapped and have total spin one. In the following we 
will consider triplets on the rungs as local triplons. 

It is natural to represent the rung states in terms of 
bond operators [30H32] , Within this representation, the 
four-dimensional Hilbert space of each rung i of the lad- 
der is given by the singlet 



4lo)<=w*=^(it;)-i;t» 



and the three triplet states 



tlx.jo), = = ^| (|tt> - 



t{i\0)i =l*o)< 

tLMi = \t+i)i 



V2 
1 



(ltt) + lW» 

(in> + ut». 



(4) 

(5a) 
(5b) 
(5c) 



The state |0) = Yl i |0),; represents the completely empty 
system. This is actually an unphysical state which is 
needed here only for formal reasons. 
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The physically most important state is the ground 
state at {x = 0, x cyc — 0} which will be used as a ref- 
erence state |ref) for our calculation. It is given by the 
product of singlets on all rungs 



is given by 



|rcf) := []l 



(6) 



One can also interpret |ref) as a singlet condensate, i.e. 
singlet operators can be replaced by unity s = 1 = 
Physically, |ref) is the triplon vacuum, i.e., the state with- 
out any excitation. 

The operator V i with a G {0, ±1} creates an elemen- 
tary excitation with S z — a (flavor) when applied to a 
singlet state on rung i. Each rung is either occupied by a 
singlet or by a triplon, resulting in the so called hardcore 
constraint on all rungs i 



sj s ■ 



€{-1,0,1}, 



(7) 



which implies particle-number conservation in terms of 
singlets and triplons. 



The Hamiltonian ( lb ) can be reformulated in terms of 
the above bond operators. One finds 



H = H ± +T + T_ 2 



To 



(8) 



where the index n in T n indicates the change of the total 
triplon number, e.g., T+2 oc tH*ss are pair creation pro- 
cesses and T_2 oc ttsts^ are pair annihilation processes. 
In contrast, To does not change the number of triplons in 
the system. This term includes triplon hopping processes 
tHs^s and two-triplon interactions tHHt. Note that we 
have omitted the spatial and flavor index here in order to 
keep the notation light. Finally, the operator H± counts 
the total number of triplons in the system. 

The low-energy spectrum of spin ladders is well studied 
[13l [HI I33H40] . All excitations can be classified by the 
parity with respect to reflections about the centerline of 
the spin ladder: an odd number of triplons has odd parity 
and an even number of triplons has even parity. 

The most relevant low-energy excitation has odd par- 
ity corresponding to a single triplon. It is completely 
characterized by the one-triplon dispersion 



uj{k) = vq + 2 v m cos(mfc), 



(9) 



m— 1 



where k is the momentum along the ladder and the v m are 
the one-triplon hopping amplitudes. The dispersion takes 
a minimum at k = tt defining the one-triplon gap A. In 
this work we are solely interested in one-triplon energies. 
Therefore, we do not discuss channels with more triplons. 
But we note that there are interesting interaction effects 
present in two-leg ladders leading to two-triplon bound 
and anti-bound states as well as pronounced continua in 
their response functions US US [36l ES I40H43] . 

The dynamic structure factor S^.'^ (k, to) is the quantity 
of interest for inelastic neutron scattering experiments. It 



S^'^(k, uj) 



— Imy^t^") , (10) 

1 — e t 



where 



X f{k,w)=i / e iwt X;e- ttr <[5?(t),^(0)])(9(*) 



(11) 

denotes the retarded dynamic susceptibility and the su- 
perscripts a,j3 € {±1,0} correspond to the components 
of the spin. Since our model is SU(2) invariant, the only 
finite components of the dynamic structure factor are the 
diagonal ones {a = (3). We define 



Srik^uj) — S^°(k 1 uj) 



(12) 



The dynamic structure factor at zero temperature is 
dominated by the one-triplon contribution [T3] . Theo- 
retically, the general form of the one-triplon contribution 
(or more generically of any single mode (SM) approxima- 
tion) reads 



S%%(k,u) = a z {k)8{u - u(k)), 



(13) 



where the prefactor a 2 (k) defines the one-triplon spectral 
weight which is known to be strongly peaked at the gap 
momentum k = tt [5J 144] . 



III. APPROACH 

In this work we are interested in the effects of thermal 
fluctuations on the one-triplon contribution to the dy- 
namic structure factor of the spin ladder. Physically, one 
expects that the triplon acquires a finite life-time upon 
heating. Then the ^-functions at zero temperature are 
replaced by resonances which are found to show asym- 
metric line-shapes [2J. 

Our approach is complementary. We do not describe 
the decay of triplons due to scattering with other ther- 
mally excited triplons, but we develop an effective single- 
mode approximation. That means that the one-triplon 
contribution to the dynamic structure factor is still given 
by a (^-function, but with temperature-dependent spec- 
tral weight. On the level of the triplon dynamics, this is 
less sophisticated than the approach by Essler and collab- 
orators [TH3]. In contrast, we are able to efficiently incor- 
porate vertex corrections which are also relevant at finite 
temperatures. This leads to a temperature-dependent 
one-triplon spectral weight in our theory. 

In the following we give a very short introduction to 
the method of continuous unitary transformations focus- 
ing on the perturbative realization which is used in this 
work. The acronym CUT stands for the continuous uni- 
tary transformations in general and PCUT for their per- 
turbative realization in particular. 
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We concentrate on the vertex corrections which are 
relevant for the single-mode description of the finite- 
temperature dynamic structure factor. For a detailed 
introduction to CUTs we refer to the literature, see be- 
low. In the sequel, we describe the mean-field decou- 
pling of the vertex corrections in detail. They give access 
to the effective temperature-dependent description of the 
dynamic structure factor at low energies. 



A. (P)CUTs 

The CUT method was invented independently by Weg- 
ner [S] and Glazek and Wilson 03 gg] in 1994. The CUT 
itself is defined through the flow equation 



d l H(l) = [r 1 (l),H(l)] 



(14) 



with the continuous parameter I € [0, oo] and the ini- 
tial condition H(l = 0) = H where H is the given 
Hamiltonian. The final form of the effective Hamilto- 
nian H c r = H(l = oo) depends on the choice of the 
infinitesimal generator ?/. 

Here we use (quasi) particle conserving CUTs [TUl EH 
l47l [48] which result in the effective Hamiltonian H c g be- 
ing block-diagonal in the number of elementary excita- 
tions q. In the case of the spin ladder studied, q is the 
number of excited triplons. In an cigenbasis \qi) (i is 
an additional index necessary due to degeneracies) of the 
particle- number operator Q with eigenvalues q € No, the 
generator ?7 pc is given by 



sgnfe - qi)Hij{l). 



(15) 



For the two-leg ladder we choose J±Q = H± giving 
an effective Hamiltonian which conserves the number of 
triplons, i.e., [H±,H c s] — 0. Consequently, the effective 
Hamiltonian H c r decomposes into a sum of irreducible 
n-triplon operators H n 



H, 



err 



oo 
n=0 



(16) 



and each n-triplon block can be treated separately, see 
Refs.HHlH] for details. 

We apply the CUTs in their perturbative version, i.e., 
the effective Hamiltonian for the two-leg ladder is ob- 
tained as a high-order series expansion in the small pa- 
rameters 



^cyc 



:= J||/Jj_ 

• — J eye I J A_ • 



(17a) 
(17b) 



The essential one-triplon dispersion cj(fc) has already 
been determined earlier |6l [12] ■ The one-triplon hop- 
ping amplitudes v m have been calculated up to order 
11 in both perturbative parameters. Various extrapola- 
tion schemes can be applied successfully [35]. Especially 



the low-energy part of the dispersion including the one- 
triplon gap can be determined very reliably up to realistic 
values x « 1.25 — 1.5 and a; cyc ~ 0.2 — 0.25. 

In order to calculate spectral properties of the system 
such as the dynamic structure factor ^(fcjw) one has to 
transform observables O by the same CUT as the Hamil- 
tonian. Hence the same generator rf c given in Eq. ( 15 ) 
is used 



d l O{l) = [rf c {l),0{l)] 



(18) 



The relevant observable 0(r) with odd parity for the dy- 
namic structure factor [5] is given locally on rung r by 

0(r) — S* t i — S^ 2 = 4,r S 0,r + to,r S 0,r- (15) 

The bare observable creates or destroys a single triplon 
with flavor on rung r. After the CUT, the effective 
observable O eS (r) comprises many channels reflecting the 
complicated nature of the interacting triplon problem. 
One finds 

^{r) = U^O{r)U = Y,°tmir) (20) 



where Of^ m (r) stands for a process where m triplons are 
annihilated and n are created. Thus n, m € Nq with the 
constraint that the difference n — m is an odd number 
reflecting the odd parity of the observable. At zero tem- 
perature, no triplons are present so that only the terms 
O n fi with n odd are relevant. 

Below we need the most relevant one-triplon contri- 
bution 0\fi which is defined by the one-triplon spectral 
weights a 2 {k). The corresponding real space amplitudes 
have been calculated earlier up to order 10 [B]. Extrapo- 
lations of the low-energy part close to k = it are reliable 
up to realistic values of x and x cyc . The contributions 
O n o with n > 3 are significantly smaller in weight than 

Oi'o BM- 

At finite temperatures, also contributions O n m with 
m > matter. These contributions give rise to vertex 
corrections. We will argue below that the most relevant 
vertex correction for the finite temperature physics at 
low energies of the spin ladder is the contribution 02,1- 
It comprises processes where a second triplon is created 
assuming another triplon is already present in the sys- 
tem due to thermal fluctuations. Thus we refer to these 
processes as conditional excitation processes. 

We determined this contribution up to order 9 in the 
perturbative parameters x and x cyc . For extrapolation, 
we used the method of internal variables [SU]. Unfor- 
tunately, no reliable Pade or DlogPade resummation on 
top of this extrapolation scheme succeeded, in contrast 
to the extrapolation of the one-triplon dispersion and the 
one-triplon spectral weights. Certainly, this fact reflects 
the more complicated nature of conditional excitations. 
Nevertheless, we expect that the extrapolation still yields 
reasonable values for the physical processes up to realistic 
values of x and x cvc . 
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B. Effective single-mode approximation 

Next we set up our temperature-dependent theory 
which results in an effective single-mode approximation. 
A single-mode contribution to the dynamic structure fac- 
tor (see Eq. (13)) depends solely on the dispersion and 



on the spectral weight. The general expression reads 



Sf M (fc,^) 



1 



1 - e 



r <4(fc)5(w-wr(fc)), (21) 



where we set the Boltzmann constant fee to unity, i.e., we 
measure the temperature in units of energy henceforth. 

First, we discuss the temperature effects in the effective 
dispersion u>T{k). Second, we incorporate vertex correc- 
tions in the spectral weights a^(fc). 



The average triplon occupation can be calculated by 
the sum over all modes (4*fc) in momentum space. At 
finite temperatures one has 



tltu) = ^TrU'J,e-T ) = 



1 



z 



dF 

duJk 



(24) 



where k is the triplon momentum. The free energy F is 
defined via the partition function Z by 



F = -Tln(Z) 



(25) 



where fee is set to unity. Integration over the whole Bril- 
louin zone yields the wanted local triplon density 



(tU) (T) = 



*k 

~2T 



dk. 



(26) 



1. Effective dispersion 

The one-triplon contribution is characterized at zero 
temperature by uj{k). In the PCUT calculation the fully- 
condensed singlet state |ref) with (s^s) = 1 is used as 
the reference state. Hence no parameter s < 1 for the 
singlet occupation needs to be introduced as it was used 
previously to take the hardcore constraint into account 
[521 |33]. The CUT takes the hardcore constraint into 
account. 

At finite temperatures, however, the situation is differ- 
ent. Even after the CUT there will be a finite occupation 
of triplons which lowers the singlet occupation. In anal- 
ogy to the previous studies at zero temperature, we take 
this effect into account by introducing the mean-field pa- 
rameter s(T) defined as the condensate value (s^) = s 
with s £ [0,1]. We stress the difference to previous stud- 
ies [32l [34] where s < 1 was used to consider quantum 
fluctuations while we use this concept only to consider 
thermal fluctations. More complex expressions such as 
s^s i are reduced to s 2 which focusses on the average value 
neglecting fluctuations. 

Then the s-dependent effective one-triplon dispersion 
reads 

U) T (k, s) = u + 2s 2 (T) v m cos(mk) (22) 

m— 1 

because each non-local hopping $ a -t a i with j ^ i implies 

the transfer of a singlet sjsj. The local processes t^ jt j 
do not require such a renormalization. 

To determine the temperature-dependent mean field 
parameter s(T) it is inserted in the hardcore constraint 
Eq. |7| yielding 

l = s 2 (T) + 3(th)(T). (23) 

So the average triplon occupation (tH) (T) determines 
s(T). Due to the spin symmetry the expectation val- 
ues (4,0 {T) with a € {—1,0, 1} are independent of a 
resulting in the factor 3 in Eq. ( 23 ) . 



If the triplons were bosons without contraint (t\.t k 

would be given by the standard bosonic occupation 
(exp(u;(fe)/T) — 1) ■ Unfortunately, there is no equiva- 
lent rigorous expression for hardcore bosons. As an ap- 
proximate treatment we follow the arguments of Troyer 
et al. [51j to reweight the partition function Z yielding 



Z = (Z u ) 



N 



Z„ 



1 



2tt 



e t dk 



where N is the number of rungs. Consequently, the free 
energy becomes 



F = -NT In 1 



2tt 



e dk 



(28) 



Combining Eqs. ( |28|26|23| yields 
3 1 



s'(T) = 1 - 



2tt Z 



e ~dk. (29) 



Due to the dependence of the dispersion wy(fc,s) on s 
the above equation defines s only implicitly. It has to be 
evaluated self-consistently. 

The value of s 2 lies in the interval [0,1]. Clearly 

lim-r^o e t = as long as uJT{k,s) > C > 
which is the physically reasonable case. Hence s 2 (T = 
0) = follows. Similarly, the limit T —¥ oo is given 

by limT^oo e t t = 1 so that s 2 (T = oo) = 1/4 is 
implied. 

For completeness, we mention another effect which im- 
plies a renormalization of the effective dispersion due 
to a finite concentration of thermally excited triplons. 
Triplon-triplon interactions of the type 44*3*4 are 
known to exist. For instance, they lead to the occurrence 
of bound states [T3J [52]. They also yield some renormal- 
ization of the dispersion upon mean-field decoupling 



t 1 t^tot 



l '2'::' l (4*3)4*4 + (*2*3)*1 *4 

+ (4*4)4*3 + (4*4)4*3 + COnst -> ( 30 ) 



G 



a) 

■ ► ■ 




Figure 2: Propagation of the excitation included in the calcu- 
lation. Diagram a) shows the one triplon excitation; diagram 
b) shows the conditional triplon excitation which requires the 
existence of a triplon. 



stands for another additional exponentially small factor. 

Restricting to the two parts O^+i m with to e {0, 1}, 
the effective observable is given by 

O off (r) = rfO(r)U 

« Ofo(r) + Of>) 

= E ( a p4,r+p + h - C -) 
V 

+ E E (<^U4.+a-+p + h - c -) • ( 31 ) 

a,/3,7 i<j,p 

The ordering of operators in the above equation is not 
unique. To avoid double counting, the convention i < j 
is used. A few excitation processes appearing in 0|i i n 
low order in the small parameters are illustrated in Fig. [3] 



which is proportional to the triplon densities \t*t) (T). 
The inclusion of the above terms will renormalize the 
effective dispersion uix(k) slightly. We do not consider 
the effect of the terms in Eq. ( 30 ) here quantitatively for 
two reasons. 

First, an estimate of the quantitative effect of these 
corrections indicates that they are not very significant. 
For instance, they are of the same order as the effect of 



the inclusion of the factor s (T) in Eq. (22 1. If this term 



is omitted no qualitative changes of our results will oc- 
cur because they are dominated by the renormalization 
of matrix elements, not by the renormalization of the en- 
ergies as we will show in the following. Second, their 
systematic treatment is very tedious because a quantita- 
tive determination of all interaction elements would be 
required. 



2. Vertex corrections and effective spectral weights 

The second building block for the one-triplon contri- 
bution to the dynamic structure factor is the one-triplon 
spectral weight. Its unconditional part, active also at 
zero temperature, is captured by Of\ after the CUT. Its 
propagation is illustrated diagrammatically in panel a) of 
Fig. p Contributions 0^ i m with m > are only active 
at finite temperatures; they represent conditional excita- 
tions. Our idea is to include such contributions as vertex 
corrections in an effective manner. We identify Ofi as 
most important because it requires only a single triplon 
to be present. Hence its quantitative contribution will 
be proportional to the triplon density. In second quan- 
tization 0|*\ contains operators of the form t^tpt^. It is 
illustrated diagrammatically in panel b) of Fig. [2] 

Any term of O^+i m with to > 1 requires the existence 
of to triplons so that its quantitative contribution will be 
proportional to the m-th power of the triplon density. 
In the regime of low temperatures not larger than about 
the spin gap each additional power in the triplon density 



♦ 

o 

♦ 



Figure 3: (color online) Some conditional excitation processes 
are sketched which are of low order in x. The arrows indicate 
the rung on which the observable O-^iM ac ts, whereas the 
cross x denotes the rung which is occupied by the one-triplon 
before the conditional excitation. Black filled circles denote 
the rungs occupied by triplons after the conditional excita- 
tion process. Open circles are in their singlet state after the 
excitation; the cyan (gray) part is shown for completeness. 

Symmetries of the system imply further constraints. 
The conservation of the total S z component for processes 
of the type ^ttt 7 lead to 



= a + /3- 7 a,/3, 7 e{0,±l}, 



(32) 



where it is understood that the bare observable is given 
by Eq. pi 



For scattering experiments we are interested in the 
Fourier transform of O (r) in Eq. (31). For the un- 
conditional contribution 0°^ we find 



Ofo(fc) 



1 J2^ r Otf(r) 

r 

— V ( e ikr n f 1 + e ikr n t 

r,p 



(33) 
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where we use 



x T (i-p) 



»(*) :=J2 e ~ ikPa P- 



(34) 



Note that the coefficients a p are real. For the symmetric 
spin ladder studied here they are also symmetric a p = 



a_ p . The operators t Q k and £o,fe are defined by 
TT^E^V 

v r 



0,k -~ 

t ._ 
O.fc -~ 



(35a) 
(35b) 



The complete treatment of the conditional excitation 
is complicated since it involves several sums and the re- 
sulting terms are not diagonal in momentum space. We 
do not, however, aim at the exhaustive description of the 
multi-particle response, but at the renormalization of the 
single-mode response. Therefore, we use a mean-field de- 
coupling for OflC^) *° identify the processes belonging 
to the creation of a single triplon. 

The mean-field decoupling is motivated by Wick's the- 
orem |53j , but represents an approximation in the present 
context because the triplons are hardcore bosons. For 
concentrations of the triplons tending to zero, i.e., for 
T — > 0, the hardcore constraint becomes ineffective and 
the use of Wick's theorem is justified 



i f J t 



(36) 



Due to the symmetry constraint ( 32 ) the expectation val- 



ues are non-zero only for the cases 



(f Jj* 7 ,p) + if = 7 and a = 
(*a,i* 7)P ) ^0 if a = 7 and = 0. 



(37) 



In all other cases the expectation values are zero. We 
define the hopping expectation values by 



r T {j-p) ■= (tl d t PtP ). 



(38) 



Due to spin rotation symmetry and translation symmetry 
this definition is independent of (3, but depends only on 
the relative distance of the two rungs j and p as depicted 
in Fig. [4j Consequently, we perform an index shift (J — 
p) — > j yielding Vr{j) which is used in the following. 



Employing the mean-field decoupling ( 36 1 renders the 
operator structure of 0|^ the same as for Of Q . There- 
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Figure 4: Illustration of the hopping expectation value Tr(j - 
p) as defined in Eq. ( 38 1 . 



r+p 



fore, Fourier transformation makes it concise 

r,i,j,p Q,/3,7 

+ e ta,r+it/3,r+jt\,r+p^ (39) 
r,t,j,p f3 

+e ikr t , r+l T T (j ~ P)<W) 

+ 7^ E ( elfcr 4,r +J T T(*-:P)^, 7 

+e ikr ^ r+3 T T (z-p)d a ^). (40) 

Shifting the indices r — > r — i and j — > j + p in the first 
sum as well as r — > r — j and z — > i +p in the second sum 
leads to 



/Jj E '^2 a if+P,p{ e lk%elkrt 0,r 



N 

r,r,3,p p 
—iki Akr 



+e-V r t 0jr ) r T (j) 



(41) 



(42) 



N 

+e- ik 3e ikr t 0ir )T T (i) 

i,j,P 

Finally the full effective observable O c (k) in momen- 
tum space is given by 

O cS (k) = Ot s (k) + Otl(k) 

= nk) + A T (k)}(tl_ k +t ok ), (43) 



where we used the successive definitions 
A T (k) := J2Mk)r T (i) 

i 

Mk) : J2 A ^ e ~' lki 



a \ ^ / 0,a,ct , a,0,a 



(44a) 
(44b) 
(44c) 



Clearly, (43) shows that the conditional excitations imply 



a renormalization of the spectral weight of the triplons 
which are the single-mode excitations for the spin lad- 
der. The temperature dependence of Ax(k) stems from 
the temperature dependence of the hopping expectation 
values Ty(j'). Similar to the calculation of the triplon oc- 
cupation in (j26j), the hopping expectation values result 
from the suitably weighted averages over the Brillouin 



mean-field (MF) approach. We demonstrate that the in- 
tensity decrease in the dynamic structure factor mainly 
arises from the contribution of the conditional excita- 
tions. The basic energy unit in this theory section is 
chosen to be the rung coupling constant Jj_ allowing for 
direct comparison for different parameter sets {x,x cyc }. 

In Ref.l2"71the parameters x — 1.5 and x cyc = 0.2 — 0.25 
were determined to describe a generic cuprate spin lad- 
der best. Thus we chose these parameters to com- 
pute the dependence of the self-consistent MF disper- 
sions shown in Fig. [5] with temperature T and wave vec- 
tor k. This provides information on how the shape of 
the effective dispersion changes upon increasing temper- 
ature and momentum. The dispersion becomes flatter 
on increasing temperature implying a larger energy gap 
A(T)=w(7r,T). 



t tU) = 



1 

2tvK 



cos(jk)dk. 



(45) 



The expressions ( 44 1 for the effective spectral weights 



hold for all momenta k. But the main effect in the spec- 
tral weight is measured at the lowest excitation energy 
A := w(7r), i.e., at the momentum k = ir. Hence we re- 
strict ourselves to k — ir in the following. For Ai{tt) we 
obtain explicitly 



Mk) = 2^ _1 < > \ a j,i+p,p + a i+pj,p 
pj 

(4}?p,p +a t+p,i,p))> ( 46 ) 



where the factor 2 comes from the fact that a ' "' 



0,-1,-1 j 1,0,1 

a„ „ and a 



-1,0,-1 i- l 

a„ ^11 „ , respectively. 



j,i+p,p 



j,i+P:P j,i+P,P , j,i+P,P ' 

Adding the unconditional and the conditional exci- 
tation processes yields the total effective temperature- 
dependent one-triplon spectral weight at k = n 



4(tt)= [a(Tr) + A T (ir)] 2 s 2 (T). 



(47) 



The factor s 2 (T) stands for the reduced availability of 
rungs for the creation of triplons if there are already some 
thermally excited. Putting all pieces together we gener- 
alize ( 13 ) to finite temperatures 



S} M (k,uj) =a^(k)S{u;-u)T(k)), 



(48) 



which we use in the present study for k — n only. This is 
the dynamic structure factor in single-mode approxima- 
tion of the dynamic structure factor including the vertex 
corrections relevant at finite temperatures. The physical 
content of Eq. ( 48 ) will be exploited in the next section. 



IV. THEORETICAL RESULTS 

Next we present our theoretical results obtained by ex- 
tending the PCUT result to finite temperature using the 
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Figure 5: Effective mean-field dispersion for x c 



0.2 and 



x = 1.5 in dependence on temperature and momentum. 

Figure [6] shows the gap as a function of temperature 
for various x = {0.75, 1, 1.25, 1.5} and a different x cyc in 
each panel. The blue (gray) dots in each panel indicate 
a characteristic temperature T c har above which the spin 
gap shows a significant dependency on temperature in 
form of a steep gap function A(T) for all x and x cyc . 
The value of T c h ar scales with the spin gap A(T = 0); 
hence it decreases with increasing x cyc . 

Three further striking points are to be mentioned. 
First, with increasing x cyc the curves for x = 
{0.75, 1, 1.25, 1.5} approach one another until they lie al- 
most on top of one another for x cyc = 0.15 — 0.2. For 
a^cyc = 0.25 they start to spread again. Second, with 
increasing x cyc the sequence of the curves in each panel 
changes from 0.75, 1, 1.25, 1.5 to 1.5, 1.25, 1, 0.75 from top 
to bottom. Third, the spin gap A(T) decreases on in- 
creasing a; C y C . 

Furthermore, it is interesting how much the average 
singlet state occupation number s 2 (T) changes with tem- 
perature. Its behaviour is illustrated in Fig. [7] Below the 
characteristic temperature T c har the occupation number 
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Figure 6: (color online) Excitation gap vs. temperature T 
for x — {0.75, 1, 1.25, 1.5} and a different x cyc in each panel. 
The red (gray) dots indicate the position of the characteristic 
temperature T cha r. 
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Figure 7: (color online) The dependence of the Average singlet 
state occupation number s 2 (T) with temperature for various 
x = {0.75, 1, 1.25, 1.5} with a different x cyc in each panel. 



s 2 (T) is almost independent of temperature for all x cyc 
resulting from the constant gap energy A(T). Clearly, 
Tchar scales with the energy gap A(0). Above T c har the 
occupation number s 2 (T) falls off the steeper the smaller 
x for a given x cyc . 
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Figure 8: (color online) The upper panel shows the MF triplon 
hopping expectation value r(j) for j — {0, . . . , 6} with x = 
1.5 and x cyc — 0.2. The lower panel depicts the associated 
conditional MF spectral weight At(it). 

The generic behaviour of the MF triplon hopping ex- 
cpectation values r(j) is depicted in Fig. [8] (upper panel) 
for j — {0, • • ■ ,6} vs. temperature for x = 1.5 and 
^cyc = 0.2. For all other pairs of x and x cyc the hop- 
ping expectation values look very similar. The associ- 
ated conditional MF spectral weight At(it) is shown in 
the bottom panel of Fig. [8] Thus, the shape of At(tt) is 
dominated by t(J) with j > 1 explaining the minimum 
in At(it) at about T = 0.514Jl which is about twice the 
spin gap A(0) [55] , 

The zero temperature one-triplon spectral weight 
a 2 (n) vs. x for x cyc = {0, 0.05, 0.1, 0.15, 0.2, 025} is shown 
in Fig. [9] With increasing x the one-triplon spectral 
weight at k = ir grows. This increasing weight on increas- 
ing x confirms the conclusion, that the most important 
contribution to the one triplon weight and to the condi- 
tional triplon excitation weight is found at k — tt. This 
is the dominant feature in the dynamic structure factor. 

Fig. [lO] depicts the conditional MF triplon spectral 
weight At (n) obtained for various x = {0.75, 1, 1.25, 1.5} 
and different values x cyc vs. temperature. The condi- 
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Figure 9: (color online) One-triplon spectral weight at k — it, 
where the spin gap A(0) occurs, for T = vs. x for various 
x cyx . The one-triplon spectral weight at k — n increases on 
increasing x for fixed x cyc . For fixed x the one-triplon spectral 
weight increases also on increasing x cyc . 
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Figure 10: (color online) Conditional MF spectral weight vs. 
temperature T for x = {0.75, 1, 1.25, 1.5} and different x cyc 
in each panel. 



tional MF triplon weight for various x shows the expected 
almost constant behaviour below the characteristic tem- 
perature T c har where all curves lie indistinguishably on 
top of one another. Above T c h ar a steep fall-off is found. 
The steepness is affected by the values of x and x cyc . The 
fall-off is the larger the larger x and x cyc are. For fixed 
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Figure 11: (color online) The momentum integrated dynamic 
structure factor I(T) = / S(k, T)dk for x = {0.75, 1, 1.2, 1.5} 
with a different x cyc in each panel is presented. It is plotted 
vs. temperature, normalized to the momentum integrated dy- 
namic structure factor Io = J S(k, 0)dk at zero temperature. 



x cyc the fall-off becomes stronger on increasing x. 

The calculated conditional MF amplitude At(tt) dis- 
plays a minimum at about J±/2 which corresponds to 
about twice the spin gap A(0). We do not see any physi- 
cal reason why the weight should increase again. Further- 
more, the justification for the approximations involved is 
less sound beyond temperatures of about 2A(0). Life 
time effects of the modes are expected to become more 
and more important [IH2] so that we do not expect the 
calculated data to be reliable for temperatures beyond 
J±/2. 

The negative conditional MF amplitude At(tt) dimin- 
ishes the one-triplon spectral weight ar(k) 2 according to 
Eq. (47) yielding the normalized mome ntu m- integrated 

I(T) I I 

structure factor -j^ 1 presented in Fig. 11 This is the 
experimentally relevant quantity. The shape of the nor- 
malized momentum-integrated structure factor in 
Fig. [TT] is very similar to the one of the conditional MF 
spectral weight in Fig. [l0| confirming that the conditional 
excitation process is the dominating effect on increasing 
temperature. This fact demonstrates the importance of 
vertex corrections in strongly correlated systems. 
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V. COMPARISON TO EXPERIMENT 

In an inelastic neutron scattering experiment the dy- 
namic structure factor is measured by the intensity of the 
scattered neutrons. In order to compare with theory a se- 
ries of measurements were made on a cuprate spin ladder 
at different temperatures T. These measurements reveal 
the temperature development of the structure factor. 

So far we discussed the theoretical results relative to 
the energy scale given by the coupling constant J±. In 
order to compare to experimental data, we have to de- 
termine this energy scale. For each pair {x,x cyc } the 
coupling constant J± is determined from the experimen- 
tal spin gap A(T = 0) = 27.6meV [23 EJ. The obtained 
values of J± are shown in Fig. 12 as functions of x (left 
panel) and of x, 

[12 



„ cyc (right panel). 
A striking point is found in the right panel of Fig. 



where for x cyc = 0.15 the value of Jj_ is nearly the same 
for all x. The sequence of the curves for various values 
of x changes at x cyc = 0.15, i.e., for lower x cyc the curves 
rise on rising x while they fall for larger x cyc . Thus the 
curve for x c 
flat. 



0.15 in the left panel of Fig. 12 is almost 
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Figure 12: (color online) Rung coupling constants J± vs. x 
(left panel) and vs. x cyc (right panel) determined from the 
experimental energy gap A(0) = 27.6meV |27l 154] . 



A. Sample and INS Experiments 

Single crystals of La4SrioCu2404i were grown using 
the 'traveling solvent floating zone' method [55] at 9 bar 
oxygen pressure. INS measurements were performed on 
the undoped ladder, La4Sri Cu24O4i using the MAPS 
time-of- flight spectrometer at ISIS, Rutherford Apple- 
ton Laboratory, U.K. The sample consisted of three co- 
aligned single crystals with a total mass of 23g (see Ref. 
[2"71 for details) . The crystals were mounted in a closed 




30 32 34 36 

Energy (meV) 

Figure 13: (color online) Background subtracted ladder sig- 
nals measured at various temperatures extracted from INS 
measurements on La4SrioCu2404i. The lower panel shows a 
zoom of the top panel. 



cycle cryostat with the (Okl) reciprocal lattice plane hor- 
izontal and the c axis perpendicular to the incident neu- 
tron beam fcj. A Fermi chopper was used to select an 
incident neutron energy of lOOmeV and 50meV and was 
rotated at a speed of 300 Hz to give an energy resolu- 
tion at the elastic line of 4.1meV and 1.6meV, respec- 
tively. Close to the gap, at an energy transfer of 30meV 
the energy resolution is 3.06meV and 0.93meV at an inci- 
dent neutron energy of lOOmeV and 50mcV, respectively. 
Data were collected at temperatures of 15K, 50K, 100K, 
and 150K, and of 15K, 50K, and 150K with an incident 
energy of lOOmeV and 50meV, respectively. Incoherent 
nuclear scattering from a vanadium standard was used 
to normalize the magnetic cross-section. In the follow- 
ing text wave vectors Q c and Q a represent the direction 
along the ladder and along the rung respectively. 



B. Background Subtraction and Data Analysis 

The ladder signal was extracted by taking a constant- 
Q cut at Q c — 0.5 and Q a — 1.2 (wave vector ranges 
0.45 < Q c < 0.55 and 0.9 < Q a < 1.5) where Q c and 
Q a are given in units of 27r/c and 2ir/a, respectively. 
The wave vector Q c — 0.5 corresponds to the position of 
the gap and the Q a range is chosen to maximise the ra- 
tio of magnetic intensity to background due to the rung 
modulation which goes as (1 — cos(Q a (i run g)), where 
rfrung is the rung distance. The background was deter- 
mined from a constant-Q cut for the wave vector range 
-0.2 < Q c < 0.2, and 0.9 < Q a < 1.5. This range 
is appropriate for the determination of the background, 
because it contains no magnetic scattering due to the lad- 
der structure factor. Each temperature run was treated 
in exactly the same way, yielding the data shown in Fig. 
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The highest intensity was measured for the coldest run 
at 15K. The experiment does not reveal any temperature 
dependence between 15K and 50K. A slight decrease in 
intensity is observed for 100K and a significant decrease 
for 150K. The spin excitation gap is found to be at E& = 
(30 ± 0.5)meV, see lower panel of Fig. [l3j which lies a 
little lower than that found in Refs. 1561 1571 and a little 
higher than found in Ref. [S7J Here we define the gap 
as the energy below the intensity maximum where the 
scattering strength is half of its maximum value. 

The FWHM of the observed peak was extracted by 
fitting a gaussian to the data. As the real shape of the 
measured curve deviates from the ideal gaussian, the fit is 
adjusted such that the gaussian and the data points with 
an energy below the energy of the maximum intensity 
match each another. The FWHM values are found to be 
5.3±0.5meV for 15K, 5.1±0.5meV for 50K, 5.3±0.5meV 
for 100K and (8 ± 0.5)meV for 150K for an incident en- 
ergy of lOOmeV and 3.3 ± 0.5meV for 15K, 3.5 ± 0.5meV 
for 50K and (5.8 ± 0.5)meV for 150K for an incident en- 
ergy of 50meV. Different FWHM values obtained for data 
collected at the same temperature but different incident 
energies was due to changes in the energy resolution. 

A broadening effect due to damping could therefore 
only be observed at 150K. This result confirms that the 
broadening with increasing temperature is small in com- 
parison to the intensity decrease with increasing temper- 
ature, which in turn validates the MF approach presented 
here. Strikingly the experimental intensity curves in Fig. 
[T3| lie exactly on top of each other above 40meV for all 
measured temperatures, indicating that there is a tem- 
perature independent tail which we attribute to the wave 
vector resolution. 

The data is compared to the temperature dependent 
theoretical spectral weight described previously. To de- 
termine the relative changes we exclusively used the in- 
tensity signal in the range 28 < lj < 40meV and cal- 
culated the integrated intensities normalized by the area 
obtained from the 15K run. Fig. |14| shows the comparison 
between theory and experiment for the temperature de- 
pendent changes of the normalized integrated intensities. 
The temperature scale is restricted to the temperature 
region that has been probed experimentally. 

Possible pairs of (x, x cyc ) matching the INS data are 
{(0.75, 0.1) , (1, 0.15) , (1.25, 0.15) , (1.5, 0.2)}. However, 
previous INS measurments |27j and optical conductivity 
measurements f6,!28! showed, that values of./- ^ 1.25—1.5 
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Figure 14: (color online) Temperature dependence of normal- 
ized dynamic structure factor ^j^j at k — -k for various pairs 
of (x, x cyc ). The energy scale is set by J± given in Tab.[I] The 
black and red (gray) dots with error bars depict the exper- 
imental result obtained from INS for lOOmeV and 50meV, 
respectively. 



^cyc) 


J ± [meV] 


(0.75,0) 


54.2 


(1.5,0) 


66.8 


(0.75,0.1) 


76.4 


(1,0.15) 


98.4 


(1.25,0.15) 


99.1 


(1.25,0.2) 


118 


(1.5,0.2) 


115 


(0.75,0.25) 


166.8 


(1.5,0.25) 


135 



Table I: Values of the coupling constant J± based on a spin 
gap of A = 30meV for the various pairs of (x, Xc yc ) for which 
data is presented in Fig. [14] 



A = 30mcV (present work) and A = 27.6meV [271154] . 
slightly different results are implied as presented in Fig. 
[Po] Both a gap of 30meV and of 27.6meV match the data 
for realistic parameters {a;,a; C y C }. In fact, using the spin 
gap 27.6meV consistent with the data from Ref. 27 leads 
0.20 — 0.27 are realistic. To be consistent to convincing agreement with the values for {x, x cyc } 

found in Ref. [37] to fit the INS data. Hence our findings 
in this work are fully consistent with previous results. 

At about half the spin gap energy A(0) our MF theory 
of vertex corrections predicts a decrease of the spectral 
weight by 25% for x = 1.5 and x cyc = 0.2 with J± — 
115meV which agrees very nicely with experiment, see 
We stress that a decrease of 25% is well in the 



and x cy 

with previous experimental results we argue that the best 
agreement within the experimental error bars is found for 
x = 1.5 and x cyc — 0.2 (black solid curve). 

In Ref.HTJcc = 1.5 and x cyc = 0.25 with J x = 124meV 
were fitted to the one triplon dispersion extracted from 
experiment. This little discrepancy in x cyc does not re- 
ally matter because of the experimental errors and the 
approximation necessary for computing At(tt). 

Furthermore, if slightly different energy gaps are as- 
sumed, such as A = 35meV (found in Refs. [55] [57]) . 



Fig. 15 



regime where we consider our approach valid and reliable. 

Note, however, that we omitted the line broadening 
due to thermal fluctuations in the present analysis which 
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Temperature [K] 

Figure 15: (color online) Normalized spectral weight based 
on different spin gaps (A = 35meV (cyan (light gray) lines), 
A = 30meV (dark lines), and A = 27.6meV (blue (dark gray) 
lines)) is shown. Dashed lines rely on (a;,a; cyc ) = (1.5,0.25) 
and solid lines on (x,x cyc ) — (1.5,0.2). The black and red 
(dark gray) dots with error bars depict the experimental result 
obtained from INS for lOOmeV and 50meV, respectively. 

play a more and more important role on increasing tem- 
perature. In particular, the unique experimental deter- 
mination of the weight becomes increasingly difficult. 



in the framework of continuous unitary transformations 
which lead to an effective description of the system in 
terms of hardcore triplons as elementary excitations. 

The vertex corrections are evaluated on the mean-field 
level, i.e., on the level of a single-model approximation. 
We found that they induce conditional excitation process 
which reduce the unconditional, zero-temperature exci- 
tation amplitude. Thereby, the spectral weight of the 
low-lying excitations is diminished because its tempera- 
ture dependence is indeed found to be dominated by the 
vertex corrections. The vertex corrections stem eventu- 
ally from the hardcore character of the triplons. Yet we 
emphasize that they go beyond the obvious mechanism 
that no triplon can be excited if its site is already excited 
by another triplon. 

We compared the calculated spectral weight quanti- 
tatively with the one measured in the undoped ladder 
La4Sri Cu24O4i using inelastic neutron scattering for a 
range of temperatures up to 150K, equivalent to half the 
spin gap energy. The theory should still apply to consid- 
erably higher temperatures as long as these stay below 
twice the spin gap energy. Experiments at higher tem- 
peratures for undoped ladders could test this framework 
further. 

Further research, both theoretical and experimental, 
on the vertex corrections in other low-dimensional quan- 
tum antiferromagnets is called for. In particular, a chal- 
lenging project consists in the combination of approaches 
studying the combined effect of line broadening of the 
single modes [IH3] and of the vertex corrections. 



VI. CONCLUSIONS 

In summary, we derived in this paper vertex correc- 
tions in the dynamic structure factor for a quantum an- 
tiferromagnet without long range order. The model sys- 
tem investigated is the two-leg spin ladder. The vertex 
corrections are linked to thermal fluctuations; they are 
relevant only at finite temperatures. We computed them 
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